arXiv:l507.08248v2 [cond-mat.quant-gas] 29 Sep 2015 


Tunable Frohlich Polarons of slow-light polaritons 
in a two-dimensional Bose-Einstein condensate 

Grusdt, F. 1 ’ 2,3 and Fleischhauer, M. 1 ’^ 

1 Department of Physics and Research Center OPTIMAS, University of Kaiserslautern, Germany 
2 Graduate School Materials Science in Mainz, Gottlieb-Daimler-Strasse fl, 67663 Kaiserslautern, Germany 
3 Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA 

(Dated: October 1, 2015) 

When an impurity interacts with a bath of phonons it forms a polaron. For increasing interaction 
strengths the mass of the polaron increases and it can become self-trapped. For impurity atoms inside 
an atomic Bose-Einstein condensate (BEC) the nature of this transition is subject of debate. While 
Feynman’s variational approach predicts a sharp transition for light impurities, renormalization 
group studies always predict an extended intermediate-coupling region characterized by large phonon 
correlations. To investigate this intricate regime we suggest a versatile experimental setup that 
allows to tune both the mass of the impurity and its interactions with the BEC. The impurity is 
realized as a dark-state polariton (DSP) inside a quasi two-dimensional BEC. We show that its 
interactions with the Bogoliubov phonons lead to photonic polarons, described by the Bogoliubov- 
Frohlich Hamiltonian, and make theoretical predictions using an extension of a recently introduced 
renormalization group approach to Frohlich polarons. 


When a mobile impurity interacts with an atomic 
Bose-Einstein condensate (BEC) it forms a polaron [T]- 
3 . These quasiparticles were first introduced by Lan¬ 
dau and Pekar 00 when they studied the electron- 
phonon interaction in polarizable crystals on the basis 
of the Frohlich Hamiltonian. One of the key predictions 
was the possibility of self-trapping of the impurity in its 
surrounding phonon cloud. The Frohlich Hamiltonian 
also provides a good description of an impurity interact¬ 
ing with a condensate, when phonon-phonon scattering is 
negligible. Using Feynman’s variational approach to the 
Frohlich Hamiltonian [6], it was predicted more recently 
that self-trapping can also take place for impurities in a 
BEC [3]. However the nature of the self-trapping in this 
system is subject of ongoing debate U . 
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FIG. 1. Ratio of the polaron mass M p to the bare impurity 
mass M as function of the dimensionless coupling constant a 
in a quasi two-dimensional BEC for different ratios of impu¬ 
rity to host-atom mass M/m. Feynman’s approach predicts 
a sharp transition for M/m < 0.01, in contrast to predictions 
from mean-field (MF) theory and an extended renormaliza¬ 
tion group (RG) approach introduced in Ref. |8|. 


For sufficiently light impurities, Feynman’s variational 
approach predicts a sharp self-trapping transition in 


three dimensions 00. Using more sophisticated theo¬ 
retical methods it has recently been claimed that, rather 
than undergoing a sharp transition, there exists an ex¬ 
tended regime of intermediate couplings before the im¬ 
purity becomes self-trapped 0 EEU H2]. In this pecu¬ 
liar regime, phonons become correlated due to phonon- 
phonon interactions mediated by the impurity. Their 
strength is determined not only by the impurity-phonon 
coupling constant a but also by the inverse impurity mass 
M -1 . We show in figH that the same is true for a 
quasi two-dimensional BEC, where Feynman’s approach 
predicts a sharp transition for ratios of impurity to host- 
atom mass M/m less than 0.01. Renormalization group 
(RG) calculations [3j 8, 12] in contrast always predict a 
smooth cross-over. 

At present, only little is known about the polaron at in¬ 
termediate couplings. Understanding this regime, domi¬ 
nated by quantum fluctuations, is of fundamental interest 
and may lead to applications in material science. For ex¬ 
ample polaronic effects may be important in the high-T c 
cuprate superconductors m, and intermediate coupling 
physics may play a role here. 

Here we propose a versatile experimental setup for 
studying polarons in a BEC at intermediate couplings 
for small impurity masses. The impurity is realized by 
coupling the condensate to a quantized mode of the elec¬ 
tromagnetic field in a slow-light (or electromagnetically 
induced transparency, EIT [14]) configuration, see FIG{2] 
(a). Here the impurity is a dark-state polariton (DSP) 
[03 EE] with an effective mass M that can be controlled 
by the control laser. We show that this tuning knob 
can be used to study the transition all the way from 
weak, through intermediate, to strong couplings. Ab¬ 
sorption spectroscopy allows to directly measure the full 
spectral function I(uj) of the polaron, from which most 
of its characteristics can be obtained HMS]. Although 
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FIG. 2. Setup for realizing tunable Frohlich polarons of pho¬ 
tons in a BEC: A quasi two-dimensional BEC of ground state 
atoms (b) coupled to lasers in a A-scheme (a). By exciting 
the driven atoms using a probe field £, a mobile impurity (a 
long-lived DSP) can be created. Its interactions with the Bo- 
goliubov phonons lead to polaron formation. The mass of the 
impurity, as well as the polaronic coupling constant, can be 
tuned by changing the Rabi frequency U c of the conrol laser. 


we concentrate here on polarons described within the 
Frohlich model, the proposed experimental setup is also 
well suited to study the impurity-BEC interaction when 
there is a sizable condensate depletion [20]. For suffi¬ 
ciently strong repulsive interactions the impurity can ex¬ 
pel the condensate from its vicinity and become trapped 
in the resulting effective potential. As argued in [ 21 ] 
this effect is particularly pronounced if the mass ratio of 
impurity to condensate atoms becomes small, a regime 
easily accesible with our scheme. 

System: We consider ultracold atoms with two in¬ 
ternal metastable states | g) and |s). They are coupled 
by a two-photon optical transition through a short-lived 
excited state \e) (decay rate 7 ), see FIG ( 2 ] (a). When 
the two-photon detuning S is within the EIT line-width, 
the non-decaying eigenmodes of this system are DSPs 
03 US], propagating with a group velocity v g much 
smaller than the vacuum speed of light Co mm- 

We assume that the atoms form a BEC in the internal 
ground state. Although v g <C Co can become as small 
as a few meters per second, it is much larger than the 
speed of sound c of Bogoliubov excitations in the BEC 
(c is of the order of a few mm/s). To avoid emission 
of Cherenkov radiation, we thus confine the longitudinal 
motion of DSPs to a single longitudinal cavity mode with 
wavenumber ko , see FIGj2](b). To minimize interaction- 
induced losses caused by scattering into excited motional 
states of atoms, we furthermore introduce a strong lon¬ 
gitudinal confinement for the atoms, leading to a quasi 
two-dimensional ( 2 D) BEC [7]. 

Now we describe how the DSPs interact with Bogoli¬ 
ubov phonons, details are presented in the supplemen¬ 
tary. The microscopic Hamiltonian l~i contains the mat¬ 
ter fields ^(r), where /i = g, s, e denotes the internal 
states and r is the transverse coordinate. The internal 
states | g) and \e) are coupled by a quantized cavity field 
£ (r), normalized such that c*£ is a 2D number density. 
$ 2 D denotes the vacuum Rabi frequency on the | g) — \e) 
transition, which is reduced by a Franck-Condon overlap 


due to the 2D confinement of the atoms (see supplemen¬ 
tary for details). The transition between \e) and | s) is 
driven by a control field of Rabi-frequency U c . 

For two-photon resonance, the DSP is given by 

fF(r) = sin# ^ s (r) — cos 6 £(r). ( 1 ) 


Up to non-adiabatic corrections, the DSP is decoupled 
from the bright-state polariton 4>(r) = cos# V ; s(^) + 
sin# £(r) which is subject to losses. Here tan# = 
# 2 Dy^Wl^c| 5 with no = Nq/L 2 denoting the 2 D BEC 
density, L being the linear system size and Nq the num¬ 
ber of atoms in the condensate. 

We assume that atoms in internal states fi and v inter¬ 
act via contact interactions with strengths < 7 ^, tunable 
by Feshbach resonances [24]. Using Bogoliubov theory 
the elementary excitations are modeled by phonons a &. 
The atomic scattering as well as the atom-light interac¬ 
tions give rise to couplings of the DSP to Bogoliubov 
phonons. We find the corresponding Hamiltonian to be 
of Frohlich [25] type (see supplementary), which forms 
the basis of all following theoretical investigations: 


T~Ly = 
+ 


COS 2 0 


= J d 2 k^u> k a ] k a k + 

J d 2 r &(r)*(r) J d 2 k e ik r V k (a h + a f _ fc ) . ( 2 ) 


Here non-adiabatic couplings to the bright-state polari¬ 
ton 4> and the excited state i/; e were neglected, but 
they are derived in the supplementary. The first term 
in Eq.Q describes free phonons, where the Bogoliubov 
dispersion is given by Uk = ckyj 1 + & 2 £ 2 / 2 . £ = 

1 /Q 

(2 mg 2 g »n 0 ) is the healing length. The speed of sound 

reads c = (^^no/m) 1 ^ 2 . The second term in Eq.Q cor¬ 
responds to the dispersion relation of a free DSP. k is the 
cavity line width and the transverse mass M of the DSP 
is determined by 


M 1 = cos 2 # M ^ 1 + sin 2 6m 1 . (3) 

Here M p h = hko /Co is the transverse mass of cavity pho¬ 
tons. The chemical potential /i is derived in the supple¬ 
mentary. The last term in Eq. (J2| describes the impurity- 
phonon interaction 


vw k 2 e \ l/4 


( k 2 e v 
\2 + k*eJ 


5eff = Sin 2 9 g 


gs • 


Vfe = 5eff „ 

Z7T \Zt /C~C“ / 

(4) 

The Bogoliubov-Frohlich Hamiltonian © is character¬ 
ized by two dimensionless numbers [ 12 j 


7 TC Z 


A 171 

and —, 

M 


(5) 


quantifying the impurity-phonon interaction and the 
mass ratio of the bosons in the condensate and the im¬ 
purity, respectively. For realistic parameters JEHU] we 
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estimate ra/M p h ~ 10 11 . By changing 6 , the mass ratio 
m/M « cos 2 # m/Mph can be tuned over a wide range. 
Typically, the impurity is much lighter than the under¬ 
lying bosons, due to its photonic component, but in the 
ultra-slow light regime mass ratios on the order of unity 
should be accessible. 

Phase diagram: The following discussion of the phase 
diagram is based on an extension of the renormalization 
group (RG) approach to Frohlich polarons introduced in 
Refs. [3] [12]. The key idea behind the earlier RG scheme 
is to decouple fast and slow phonon degrees of freedom 
perturbatively in every momentum shell. In Ref. [8] we 
extended this approach by performing a global mean- 
field (MF) shift after every RG step, corresponding to 
an inclusion of infinitely many diagrams. The extended 
method is not only more accurate for strong couplings, 
but it is also necessary to calculate the effective polaron 
mass in a regime where the impurity is light 0. 

In FIGj3]we present the full phase diagram of the 2D 
Bogoliubov-Frohlich polaron. We distinguish three dif¬ 
ferent regimes of weak-coupling (where Lee-Low-Pines 
MF theory [27] applies), strong coupling (where Lan¬ 
dau and Pekar’s strong coupling approximation applies 
00) and intermediate coupling (where neither of the 
two approaches is accurate). All regimes are connected 
by smooth cross-overs, and we estimated the transition 
points from the behavior of the effective polaron mass as 
detailed below. It can be shown analytically (see Ref. [8]) 
that MF theory is not only asymptotically exact in the 
commonly discussed limits a 0 and M oo, but also 
in the limit where M —>> 0. 
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FIG. 3. Full phase-diagram of the two-dimension Bogoliubov- 
Frohlich model: For sufficiently light impurities and large 
enough a an extended regime of intermediate couplings is 
found. Data points were estimated from curves as those 
shown in Fig|lJ Parameters used in the RG simulations were 
Ao = 2000/£ and P = O.OlMc (for their definition see Ref. [8]). 
We also plotted the maximum values awx below which the 
Frohlich model is valid, for rj = 1,10 as defined in the text. 
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In the case of BEC polarons realized with bare atoms 
different polaron regimes can be accessed only by tuning 
a while the mass ratio is fixed around a value between 



FIG. 4. DSP-polaron mass Mp as a function of the inverse 
transverse mass M _1 of the DSP. Results from the extended 
RG scheme introduced in [8] are compared to predictions from 
MF. Parameters are Ao = 2000/£, P = O.OlMc, and a = 158. 


~ 0.1 to ^ 10. For the photonic setup, on the other hand, 
the impurity mass M can be used as a tuning parameter. 
In particular, extremely light impurities can easily be 
created and regimes of the phase diagram inaccessible to 
bare atoms can be addressed. This versatility makes the 
photonic setup superior to purely atomic systems, for 
the investigation of the transition from weak, through 
intermediate to strong couplings. 

Now we turn to a more detailed discussion how the 
phase diagram in FIG(3]was obtained from the extended 
RG approach [8]. We expect that in the suggested ex¬ 
perimental setup our calculations can be put to a test. 
In FIGjlJwe show an example how the effective polaron 
mass depends on a for a light impurity ( M/m = 0.01). 
As found previously using the perturbative RG 0HD , a 
smooth cross-over takes place from a quasi-free polaron 
to a self-trapped polaron. For small couplings, M p in¬ 
creases linearly with a according to the MF prediction 
and crosses over into the intermediate coupling regime 
with a non-linear growth. Eventually the strong-coupling 
regime is entered where M p increases linearly with a 
again, but with a different slope. We use these crite¬ 
ria to define the different regimes in the phase diagram, 
shown in FIG [3] 

In FIG.[4]the effective polaron mass is shown as a func¬ 
tion of M/m for large interaction. One recognizes a very 
substantial increase of the polaron mass for light impu¬ 
rities, much larger than predicted by MF theory. We 
expect such dramatic effects of quantum correlations to 
be easily detectable in the proposed photonic setup. In 
the limits of extremely heavy impurities (i.e. M —»• oo) 
as well as extremely light impurities (i.e. M —)> 0) on 
the other hand, the MF prediction becomes asymptoti¬ 
cally exact and the polaron mass approaches that of the 
impurity. In the first case the phonon dressing cannot 
influence the impurity mass and in the second case the 
impurity is too fast for the phonon cloud to follow. 
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Experimental considerations: DSPs in ultracold 
BECs have been observed experimentally in the slow- 
light limit HI HU- By performing similar experiments in 
quasi two-dimensional BECs [7 with light confined to 
a cavity, see FIG [2] photonic Frohlich polarons can be 
realized. By varying the intensity of the control laser U c , 
the effective mass of the DSP can be tuned, and using 
Feshbach resonances [24] the coupling strength a can be 
varied. This should allow to explore the phase diagram 
shown in FIG [3] Realistic experimental parameters are 
provided in the supplementary material. 

Next we discuss conditions when the Frohlich Hamil¬ 
tonian is valid. In the derivation of the model [3CE2] 
we neglected phonon-phonon scattering induced by the 
impurity, which is justified when e := ynph/^o <C 1. 
Here n p h denotes the (real-space) phonon density [2, |3] 
which we estimate by n p h ~ Aph/£ 2 . The phonon num¬ 
ber (at zero total momentum P = 0) can be calcu¬ 
lated from MF theory and the condition e<l becomes 
e = « 1- Demanding an upper bound 

e < e max thus constrains |g e ff | which yields an upper 
bound for the coupling strength, 


L max 





Here 4 a gg is the extend of the quasi two-dimensional 
BEC in the strongly confined region and a gg is the scat¬ 
tering length of ground state atoms (in three dimensions). 

To estimate which range of parameters in the Frohlich 
polaron phase diagram can be accessed in an experiment, 
we plotted <a max for p = (4 /a gg ) = 1 and 10 

in FIG [3] One recognizes that the validity of the Frohlich 
model for BECs extends into the intermediate coupling 
regime, while going to strong coupling may require to 
go beyond the Frohlich model. Ultimately experiments 
need to clarify how the system behaves at intermediate 
couplings, and we believe that the proposed setup is well 
suited to explore this. 

Experimental signatures: We proceed by discussing 
possible signatures of polaron formation. In order to cre¬ 
ate a DSP-polaron one can envision to first store a weak 
probe pulse, ideally containing a single or a few photons, 
in the BEC using the storage protocol of [29] 30J and 
subsequently restore an intra-cavity DSP with small pho¬ 
tonic component, i.e. 0 < cos 2 0 <C 1. Most strikingly, 
the effective mass of the polaron M p significantly in¬ 
creases as compared to the bare mass M, see FIGg] One 
way to measure this effect is to observe dipole oscillations 
m of a polaron wavepacket inside a harmonic potential 
Muj^r 2 /2 seen by the DSP. The weak harmonic confine¬ 
ment with an oscillator length i = (M p cjo) -1 ^ 2 £ can 

easily be implemented using spherical cavity mirrors. 

A more powerful method for analyzing photonic po¬ 
larons is absorption spectroscopy upon driving the cavity 
by an external laser at frequency uj and with momentum 


P, i.e. with amplitude £ ~ £ 0 e lP ' r ~ luJt . The absorption 
rate T of photons from the laser is given by the spectral 
function /(cj, P), T(o;,P) ~ /(cj,P), in complete anal¬ 
ogy to the radio-frequency spectroscopy discussed e.g. in 
Ref. pi] . 

The momentum-resolved spectral function of the pho¬ 
tonic polaron has a characteristic delta-function peak 
I c oh4b P) = Z5(uj — E 0 (P)) which is located at the po¬ 
laron energy uo = Eq{P). By measuring the momen¬ 
tum dependence of the polaron energy (around P = 0) 
the polaron mass can be obtained. Using the sum-rule 
f duj /( cj , P) = 1, also the quasiparticle residue Z can be 
obtained from the spectral function. 

Summary: In this article we suggested a realistic ex¬ 
perimental setup for exploring polaron formation of mo¬ 
bile impurities inside a BEC. By coupling the atoms to 
lasers in a slow-light setting, we showed that DSP impu¬ 
rities with a tunable mass can be realized. Their interac¬ 
tion with the Bogoliubov phonons of a BEC can be mod¬ 
eled by a Frohlich Hamiltonian. One of the main motiva¬ 
tions to study this system is to explore the self-trapping 
transition experimentally, with the impurity mass serv¬ 
ing as a flexible tuning parameter. The physics of this 
transition, dominated by phonon correlations, is poorly 
understood. The theoretical analysis presented suggests 
a smooth cross-over rather than a sharp phase-transition 
as may be expected from Feynman’s variational approach 
[3!. Ultimately, experiments are needed to clarify how the 
polaron becomes self-trapped. 

The suggested setup furthermore raises new questions, 
including how the polaron properties change in a regime 
where the Frohlich Hamiltonian is no longer valid and the 
formation of bubble polarons may be expected [ 21 ]. Al¬ 
though the validity of the theoretical analysis presented 
here is questionable in this parameter range, the ability 
to tune both the coupling strength a and the impurity 
mass M at will makes our system appealing to search for 
new physical effects. 

Finally in solid state systems polarons have almost ex¬ 
clusively been studied under equilibrium conditions. Ul¬ 
tracold quantum gases provide long coherence times and 
allow to study dynamical effects. This includes the possi¬ 
bility to measure the full spectral function [13 HU, which 
is possible in our system using absorption spectroscopy. 
Also the dynamics of polaron formation can be studied in 
real time in the suggested experiments. The use of pho¬ 
tons coupled to short-lived atomic states moreover opens 
the possibility of studying polarons in driven-dissipative 
systems far from equilibrium. 
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Supplementary Material: Tunable Frohlich Polarons of slow-light polaritons in a 

two-dimensional Bose-Einstein condensate 


DARK-STATE POLARITONS IN A BEC 

In the following we describe in detail our formalism for the description of photonic polarons in a BEC. From the 
full microscopic model we derive an effective Frohlich polaron Hamiltonian for the long-lived DSPs. 


Microscopic Hamiltonian 

The starting point for our analysis is the 3D Hamiltonian describing free bosonic atoms of mass m in their three 
internal states | g), \e) and |s), as well as the light-field, the atom-light interactions and inter-particle scattering, 

H = 7io + Hal + Hint. (7) 

In the absence of interactions between the atoms, the excitations of interest are described by 

tio = Jd 3 x^l(x) + ^( x) + + ip e (x) 

+ #(*) (“2 +<5 + K ° (z) ) ^(*) +^(*) - iK ) ^(*)}> ( 8 ) 

where we made use of a rotating frame oscillating with the frequencies of the optical fields £ and £ 2 C , respectively. 
We introduced the one and two-photon detunings A and S which determine the energies of the atomic states in the 
rotating frame, see FIG. 2 (a) in the main text. Moreover we introduced the decay rate 7 of the excited atomic state, 
and the cavity loss rate k. 

In Eq.(| 8 | we consider a single longitudinal optical mode of the cavity. The corresponding mode function for a Fabry- 
Perot resonator is £ z (z) = y/2]L z cos(koz), where L z is the length of the cavity and ko the longitudinal wavenumber 
of the mode. The dynamics of photons in this mode, 

£{x)=£{r)£°{z), x = (r,z), (9) 

is determined by the transverse part r only. The effective transverse photon mass is M p h = hko /Co with Co denoting 
the vacuum speed of light. Similarly, the atoms are trapped inside the cavity and we assume that they are confined 
to the lowest mode of a strong axial potential V®(z) (created e.g. by additional laser beams). Their annihilation 
operators can thus be written as 

$t(n)=$t(r)ipl(z), (j, = g,s,e, ( 10 ) 

where the normalized longitudinal mode functions ^(z) are determined by the ground state solution of the Schrodinger 
equation (—<9^/2 m + V®(z)) ^(z) = By integrating out the longitudinal direction 2 , we obtain an effective 

two-dimensional model, 

n 0 = Jd 2 r^l(r) (-^) $ g (r) + ft e (r) - *(7 + *A + is° e - </> e (r) 

+ 4>t(r) ^-^+<5 + e°-£°^ s (r-)+f t (r) ~ £( r )}> ( n ) 


where we chose a frame rotating with the additional frequency s®. 
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The interactions of the atomic states with the cavity light field £ and the classical control field f^ c , see FIG. 2 (a) 
in the main text, are described by 

■Hai = (12) 

We already eliminated the longitudinal degrees of freedom in this expression, as described in the previous paragraph. 
This results in an effective two-dimensional coupling constant, 

g2D=gf dz (13) 

Jo 

with g being the vacuum Rabi-frequency of the | g) — |e) transition. We assume here that the cavity mode has no 
coherent amplitude, i.e. (£) = 0. The following analysis can be generalized to the case when a finite amplitude £q is 
present corresponding to £ £ + £q. We will remark on the effects of this at the end of this section. The atomic 

interactions read 

"Hint = J d 2 r j^- (yt(r)) (y,(r)) + 4 ?^( r )V’s( r )V’l( r )^( r ) 

+ gge^l( r Hg( r )^l( r )^e(r) + g^^l(r)ip s (r)^l(r)ip e (r) j, (14) 

again after the elimination of longitudinal degrees of freedom. The effective two-dimensional interaction strengths are 
given by g j® = g^ f Q Lz dz \ip°(z)\ 2 \^°(z)\ 2 . 


Derivation of the Frohlich Hamiltonian for DSPs 


To simplify the Hamiltonian 0 we first treat only its quadratic part, coupling the matter fields Vv( r )- We start 
by considering the case when both the one and two-photon detunings So = S + — s® and Ao = A + — s® in the 
free Hamiltonian vanish, = Aq = 0. 

We assume that a macroscopic number of atoms is initially prepared in the system. They are optically pumped 
into the atomic ground state on a time scale r g = ^/\Vt c \ 2 . When being cooled down they condense in the ground 
state and form a BEC, provided the thermalization rate is sufficiently large. 

Next we describe the macroscopic condensate of Ao atoms in the ground state using standard Bogoliubov theory, 
see e.g. Ref.pQ. We can now write 

^ g (r) = ~(4, o\/No + U k a k - v k CL ] _ k j , (15) 

k 


and as a result the effective Hamiltonian for ground-state atoms becomes = ^2 k w k Q> k Q> k + E eec (No). Here a k 
annihilates a Bogoliubov phonon with momentum k and E eec (No) is the macroscopic ground state energy of the 
BEC. L denotes the linear system size in transverse direction. The Bogoliubov form factors are given by 


u k 1 _ 1 / 2m + ffDD n 0 

Vk J \f2 V U k 


(16) 


where no = Nq/L 2 is the two-dimensional density of the BEC. 

Now we replace the ground state field operator ^ g in terms of the condensate fraction y/NoS k ,o and the Bogoliubov 
phonons a k using Eq. Quantum fluctuations in the system, including the impurity, are thus described by the 
field operators a k , ^ s (k), ^(fe) and £(k). Before treating their mutual (non-linear) interactions, which give rise 
to the effective Frohlich Hamiltonian, we discuss terms in the Hamiltonian which are quadratic in these quantum 
fluctuations. We obtain 


7/fluc = 


= j d 2 k(s\k) -in) £(k)+fy(k) (^ - 27 + A 0 - /z B ec^ 4>e(k) 

+ (^ + n 0gl? - A^bec) i>s(k) + (n c fteik^sik) - J2dV«o i>l(k)£(k) + h.c.) (17) 
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where we introduced the chemical potential /tbec = E eec (Nq) — E eec (Nq — 1) = n 0 g 2E of the BEC. 

Now we solve 'Hr uc approximately by introducing dark (4/) and bright (<f>) state polaritons, 

T(r) = sin# ^ s (r) — cos 0 £(r), 4>(r) = cos# ^ s (r) + sin# £{r). (18) 

With the atom-light mixing angle being defined as 

tan# = ^dv^oM, (19) 

the effective Hamiltonian reads 


— - i'y + A 0 - p-bec ) 'ipeik) 


Hn uc = J d 2 k[¥(k)*(k)v™r + &( k )$(k)v* SP +iPt(k) 

^l(k)^(k) g| D no cos 2 d +^\k)&(k)£ly<f>(k) + h.c. J | (20) 


The dispersion relations of the dark and bright state polaritons are 


,dsp 


,bsp 


k 2 


Vv, = 


2M ph 

k 2 


2 M 


ph 


k^ 
2m 
k 2 
2m 






— Mbec^ 

) sin 2 #, 

(21) 

— Mbec J 

| cos 2 #, 

(22) 


and the non-adiabatic coupling between dark and bright polaritons is given by 

k 2 k 2 


Q\fr&(k) = cos# sin# 


+ in + n 0 g 2 ® - ubec I • 


^gs 


(23) 


v 2m 2M p h 

From Eq.(|20|) we observe that the DSP is a long-lived excitation of the system, which is only weakly coupled 
to the bright polariton via for cos# <C 1. The bright polariton in turn is strongly coupled to the short-lived 
excited state, \Q& e \ Therefore we can restrict our analysis to the study of DSPs and neglect non-adiabatic 

couplings to other states (4> and to be precise). Non-adiabatic corrections are discussed in the following section of 
this supplementary material. 


Expressing the remaining terms in the interaction Hamiltonian (14) using the polariton basis (consisting of 4/, 4>, 
i/j e and a*.), some algebra yields the effective Frohlich Hamiltonian for DSPs interacting with the Bogoliubov phonons 
in the BEC, see Eq. (2) in the main text. The effective coupling strength is given by g e q = sin 2 0g 2E . We ignored all 
scattering channels into the remaining fields (4> and i/j e ) and neglect corrections to the Frohlich Hamiltonian resulting 
from the atom-light interaction. This is justified for realistic system parameters. In addition, two-phonon scattering 
on the spin state \s) were neglected, giving rise to the condition in Eq. (6) of the main text |2] 3]. 


Cavity field with coherent amplitude 

In the above discussion we considered a DSP formed by a quantized few-photon cavity field. In general the cavity 
field may have a non-vanishing coherent amplitude £o, i.e. £ has to be replaced by £ + £q. In this case the derivation 
of the effective Frohlich Hamiltonian is a bit more involved and leads to small modifications of the coupling constants 
as well as to some additional correction terms. We do not present the full derivation here, but only remark on 
the main additional features. First of all the simultaneous presence of two coherent fields, Q c and £o, drives the 
atoms into a dark state of the internal dynamics, which is a superposition of ground-state and spin-state atoms, 
'ipD(r) = ski'd V’s(r) — cosd ^ g (r). The orthogonal bright state, described by ^b{t) = cosd ^ s (r) + sind ^ g (r) is 
strongly coupled to the short-lived excited state ^(r). Here the atomic mixing angle is determined by 

tand = 02 d|£o|/|^c|- (24) 

For sufficiently low temperatures the atomic gas forms a BEC in the dark state ^d- In this case the dark and bright 
state polaritons are generated by quantum fluctuations of the cavity field £ and the atomic bright state ^b 

T(r) = sin# ^b(^) — cos# £ (r), 4>(r) = cos# t/S b( r) + sin# £{r). 


(25) 
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As a consequence of the mixture of ground- and spin-state atoms all interaction strengths become renormalized. 

52D ->■ 52D COS 2 ■& (26) 

gf 9 -»• 9l D D = gf g cos 4 V + <?L D sin 4 0 + 2g™ sin 2 0 cos 2 0. (27) 

g™ -> 5 2 ? = 2 cos 2 sin 2 (*£ + 5 2 ? - 5 2 ?) + 5 2 ? 

and some new scattering and non-adiabatic coupling processes arize, which we do not list here in detail however. In 
the limit of a weak cavity field, |^2 d|^o he. sin 2 $ <^C 1, the same effective Frohlich polaron Hamiltonian as 

given in Eq.(l) of the main text can be drived. 


NON-ADIABATIC CORRECTIONS TO THE DSP PICTURE 


Now we discuss further under which conditions the effective Frohlich model for DSPs is valid. This part extends 
our derivation presented in the previous section. Here for simplicity we restrict ourselves to the case when £ 0 = 0 
vanishes, such that ski'd = 0. We calculate the effect of non-adiabatic mixing of dark- and bright polaritons via 
discuss corrections to the Frohlich Hamiltonian, and derive interaction-induced DSP loss rates. 

So far we treated the cavity decay ^ ft, the kinetic energies ^ k 2 and the effective detuning 5 e ff = — //bec 

perturbatively to leading order, leaving the definition of polaritons unaffected. Different methods were employed 
previously to obtain results from higher-order perturbation theory for the polariton dispersions 0E]. Here, however, 
we are not only interested in the corrections to the polariton dispersions ft^ SP , but also in the corrections to the actual 
eigenstates. These give rise to modified interactions with the Bogoliubov phonons. For example, when the dark-state 
polariton ^ acquires an admixture of the excited atomic state V’e? the non-linearities of the atom-light interaction 
give rise to additional terms in the Frohlich Hamiltonian. 

As a first step towards obtaining non-adiabatic corrections, we introduce a new basis rotation 


Xi 

X2 


1 

71 


-i i 
i i 



(28) 


To a good approximation this diagonalizes the coupling D$ e between 4> and x/j e , because under the slow-light conditions 
we are interested in it holds Q$ e ~ # 2 Da/^o |T|, where T = 7 + zAq — ig bec- The new Hamiltonian now reads 


Tiflu 


J d 2 fcj^ SP &(k)9{k) +vf x\{k)xi{k) + v% 2 xl(k)X 2 (k) 


2 TO 


* r -"fc SP ) (x\(k)x 2 (k) +xt(k)xi(k) S j 


+ 


(&) 

~VT~ 


(xUk)+x\(k))*(k)+h. c. 


(29) 


All off-diagonal couplings in this Hamiltonian are small compared to D$ e , scaling with T, ft, k 2 or S e q. The dispersions 
of the new bright polaritons are given by 


ft ? 1 = - 

k 2 


1 (tf_ 

2m 


- iT + ftj 


BSP 


- 2 D 


<£e I 5 


ft ? 2 = - 

k 2 


1 ( 

2 m 


- iT + ft; 


BSP 


2 a 


<£e 


(30) 


Now we calculate higher order non-adiabatic corrections to the dark-state polaritons \I/, caused by the couplings 
Cl^(k). In this step not only the energies but also the basis states are modified by the perturbation. To this end we 
employ a generalization of the Schrieffer-Wolff transformation to non-Hermitian systems. We define a new basis by 


xi 

X 2 


/ 1 Q _ 

v/ 20 ^ e 

Q ^ 

j 

\/20$ e \/20<j> e 



(31) 


up to corrections of order D^ 2 . 

In the new Hamiltonian, the dark-state polariton is approximately decoupled from both bright polaritons Xi^ 
(up to terms of order Dj 1 on the off-diagonal of the effective Hamiltonian), 


"Hfluc = Jd 2 k <JA SP #' t (fc)'P(fc) +V* 1 Xi ( k )x'i(k) + pf ic 2 (k)x' 2 (k)+ 

~ 2 - ir -^ SP ) (x'd^W+^WxiW)}- (32) 
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There are no corrections to the DSP dispersion, A^ SP = 0, because corrections from x[ and x' 2 cancel each other to 
leading order in Dj*. (If, however, we conduct perturbation theory in r/U$ e , we checked that the same corrections 
to the DSP dispersion are obtained as reported previously in Refs. 120 ) 

Finally we are in a position to calculate the effective DSP-DSP interactions. By introducing the new basis x[ 2 , 4/' 
we obtain to following correction to the Frohlich Hamiltonian, 


AH f = 1 ^^' ,t (fe + q)# / (q)| cos 2 9 sin 9 (? n 0 g 2 ° - Mbec + *k) 

k,q l 


k 2 /2m 


2g?9 n o + k 2 /2m 


\ 1/4 

) (®fc + + 


02D 


+ ~'Z‘ J cos 2 9 sin 9 [ —— 


2 fi. 


<E>e 


1 


m Mr 


P h 


u k 


(k + q) 2 a k + q 2 al_ k -v k (fe + q) 2 a J _ k + q 2 a k 


2~t 


(33) 


The term in the first line modifies the effective interaction strength g e ^ in the Frohlich Hamiltonian. The terms in 
the second line depend explicitly on the momentum q of the impurity, yielding terms which are not contained in the 
usual Frohlich model. 

Our formalism readily yields all interactions between the fields x' 12 , ^ and phonons as well, but the expressions 
are too cumbersome to write them out explicitly. Here, instead, we consider only the leading-order scattering terms 
causing interaction-induced DSP losses (and a renormalized DSP dispersion). Within the Born-approximation only 
the following two terms are relevant, 


AHint = E E { W k )& - k xl(Q + *)*(<*) + h.c.} 

k,q 1= 1 


(34) 


The corresponding scattering amplitudes are given by 

„2D 


w£ 1,2) = ^cosflsinflkp f —) Tv k -^-cos9, 

k y/2 L 2 \2g 2 »n 0 + k 2 /2m ) k y/2L 


(35) 


where the ” (the ”+”) sign corresponds to £ = 1 (to £ = 2). From the Born approximation we obtain the 
renormalization of the DSP dispersion, 


A^ sp = - 


2 r 

£/ 

1=1 J 


d 2 k 


\ W l 


09 12 


U k + V h 


X& _ v DSP' 


(36) 


Therefore the interaction-induced DSP losses scale like 7dsp/ 7 (<7p?) 2 c °s 2 0, (#2d) 2 cos 2 0, up to corrections of 


order cos 2 0. 


REALISTIC NUMBERS 

Finally we we want to provide some realistic numbers for all model parameters. We consider the hyperfme states 
of 87 Rb, | g) = \F m 1 = —1), |s) = | F m 2 ,mp = 1) with scattering lengths a gg = 100.4ao, a gs = 98.006ao, 

a ss = 95.44ao [6] where ao is the Bohr radius, which are coupled in a Lambda scheme with ko = 27r/780nm (D2 
line). Furthermore we assume a density of no = 100 fim~ 2 for a transverse confinement characterized by the oscillator 
length i z = 0.29/rm as in the experiment of Ref. [7 . We assume a cavity length of L z = 30/rm, estimate typical 
atom-light coupling strength g^D ~ 10 x 27r m 3 / 2 /s [ 8 ] and use a small Rabi frequency U c « 20 x 2tt 10 3 /s (the 
relation to two-dimensional quantities is discussed in the methods). As a result we find M/m ~ 2 and a ~ 0.03. By 
increasing the Rabi frequency to e.g. U c « 20 x 2i r 10 5 /s a much smaller mass ratio of M/m ~ 2 x 10 -4 can be 
realized. The coupling strength a can be tuned via Feshbach resonances, and we find that the Frohlich Hamiltonian 
is valid below a max « 1.4 (<a max ^ 5 x 10 3 ) for D c « 20 x 2tt 10 3 /s (U c ^ 20 x 27r 10 5 /s). This value is close to where 
the intermediate coupling polaron regime begins. 
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